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Abstract 

We study the induced photon bremsstrahlung from a fast quark produced in AA-collisions 
due to multiple scattering in quark-gluon plasma. For RHIC and LHC conditions the in- 
duced photon spectrum is sharply peaked at photon energy close to the initial quark 
energy. In this region the contribution of the induced radiation to the photon fragmen- 
tation function exceeds the ordinary vacuum radiation. Contrary to previous analyses 
[4, 5, 6, 7] our results show that at RHIC and LHC energies the final-state interaction 
effects in quark-gluon plasma do not suppress the direct photon production, and even 
may enhance it at px ~ 5 — 15 GeV. 

1. In resent years much attention has been attracted to the direct photon production 
in v4v4-collisions at RHIC and LHC energies (see, for example, [1, 2, 3] and references 
therein). It is expected that for sufficiently small transverse momenta (p? ^3 — 4 GeV) 
the dominating source of the direct photons at RHIC and LHC is radiation from quark- 
gluon plasma (QGP), and at higher p? it is hard partonic mechanisms (Compton-process, 
quark- ant iquark annihilation, and bremsstrahlung from fast quarks (antiquarks) produced 
in hard reactions) [1, 2]. Nuclear effects should modify the pQCD partonic contribution 
to the direct photons in AA-collisions as compared to that in pp interaction. For the 
Compton and annihilation processes which occur at small distances this modification, 
related to the initial state interaction (ISI) effects (nuclear shadowing and Cronin effect), 
is relatively small. However, one can expect strong final-state interaction (FSI) effects 
in the QGP for the bremsstrahlung which occurs at large distances. Investigation of the 
influence of the FSI on photon bremsstrahlung is of great interest from the point of view 
using the direct photons as a probe for the QGP at RHIC and LHC. It is especially impor- 
tant for LHC energies where the bremsstrahlung component is the dominating partonic 
mechanism [1, 4, 5]. 

It was suggested [4, 5, 6, 7] that at RHIC and LHC the quark energy loss in the 
QGP phase due to the induced gluon radiation related to the multiple scattering should 
suppress strongly the bremsstrahlung contribution. At LHC energies it is equivalent to 



strong suppression of the total contribution of the pQCD mechanisms. In [4, 5, 6, 7] it was 
assumed that the only effect of QGP on photon bremsstrahlung comes from the shift of the 
quark energy to a smaller value due to gluon emission before photon radiation. However, 
the analyses [4, 5, 6, 7] missed two essential points. First of all, multiple scattering which 
fast quarks undergo in the QGP must enhance photon radiation due to the induced photon 
emission. Also, the assumption that all gluons are radiated before photon emission is not 
justified since the formation lengths for gluon and photon radiation are of the same order. 
In this case the gluons radiated after the photon do not suppress photon emission. For 
this reason one can expect that analyses [4, 5, 6, 7] overestimate suppression of photon 
bremsstrahlung. In the present paper we address the photon bremsstrahlung at RHIC 
and LHC energies accounting for the induced photon radiation and the effect of finite 
gluon formation length. 

2. The contribution of the bremsstrahlung mechanism to the cross section of photon 
production can be written as [8, 9] 



^ A (Pt) [dx da AA (p T /x) 

-^D q ^(x,p T /x) q , (1) 
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where da AA {pr) / dydp\ is the cross section of the processes A + A — > q + X (hereafter 
we suppress the argument y and mean the central rapidity region y pa 0), D q ^ n (x, E) 
is the fragmentation function for the q — > jq transition which accounts for all the FSI 
effects (x = Ej/E, E 1 and E is the photon and initial quark energies, respectively). 
Summation over quark (antiquark) states is implicit on the right-hand side of (1). We 
write D q ^(x, E) in the form 

n / 77i\ _ dP va c( x iE) dPi n d(x,E) , 

E) ~ Tx + Tx • (2) 

The first term on the right-hand side of (2) is the probability distribution of the ordinary 
vacuum q — > jq splitting and the second one corresponds to the induced transition due 
to quark multiple scattering. 

An accurate evaluation of the induced photon emission from a fast quark requires 
treatment of photon and multiple gluon radiation on an even footing. We begin with 
discussion of the induced q — > jq transition ignoring gluon radiation. A qualitative 
method for accounting for the gluon effects will be discussed later. To evaluate the induced 
spectrum we use the light-cone path integral approach [10] (see also [11, 12, 13, 14]) which 
allows one to account for the Landau-Pomeranchuk-Migdal (LPM) effect [15, 16] and 
finite-size effects which play an important role in the problem of interest. The induced 
photon spectrum for a fast quark produced at z = can be written in a form similar to 
that for gluon emission given in [17] (we take z axis along the quark momentum) 

dPind f j_ i _ s d&eff (x, z) 



dx J dx 
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where n(z) is the number density of the medium, and 
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is the in-medium (z-dependent) Bethe-Heitler cross section. Here a(p) is the dipole cross 
section of a quark- ant iquark pair of size p with a particle in the medium. ty(p,x) is the 
light-cone wave function for the q — > 75 transition in vacuum, and fy m (p, x,z) is the in- 
medium light-cone wave function at the longitudinal coordinate z (we omit spin indices). 
The no spin flip wave functions, dominating the spectrum, read 
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where P(x) = ie q ^a em /2x[s y (2 — x) + 2s g a;]/2M(a;) (s ?;7 denote quark and photon helic- 
ities), Lj = 2E(1 — x)/m 2 q x is the photon formation length, JC is the Green's function for 
the two-dimensional Hamiltonian 



H(z) = - 
with M(x) = Ex(l — x), and 
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is the Green's function for the Hamiltonian (7) with n(z) = 0. 
The dipole cross section reads 



a(p) = C(p)p 2 , 



where 



C(p) 



CtC f 



exp(iqp)] 



(9) 



(10) 
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Here Ct,f is the color Casimir for the medium constituents (quarks and gluons) and 
quark, p£> is the Debye screening mass. 

The effective Bethe-Heitler cross section (4) differs from that for a quark incident 
from infinity on an isolated scattering center due to the LPM effect and finite-size effects 
originating from the in- medium light-cone wave function entering (4). In the high-energy 
limit, when L <^ Lj, here L is the thickness of the medium, the typical values of p 

become small (< 



<2L/Ex(l — x)), and the LPM effect can be neglected. In this case the 
spectrum is dominated by the N—l scattering. In this regime the most important effect 
is modification of the in-medium light-cone wave function due to the finite-size effects. 
Let us first discuss the N — l term to illustrate qualitatively the role of the finite-size 
effects. Neglecting the Q 2 -dependence of a s from (4)-(10) one can obtain 



d<J* f H f(x,z) 
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(hereafter for simplicity it is implied that e q = 1). The derivation of (11) is similar to 
that for gluon emission given in [18]. After substitution of (11) in (3) one obtains (for 
n(z) =const) the spectrum oc L 2 

dP ind (x, E) _ -Ka em a s C T C F nL 2 {l - x + x 2 /2) 

dx ~ 8E(l-x) ' [ ' 

One sees from (12) that in the high-energy limit contrary to the ordinary Bethe-Heitler 
spectrum cx 1/x the bremsstrahlung is oc 1/(1 — x). Of course, formulas (11), (12) 
become invalid at (1 — x) <C Lm 2 /E when <^ L and the spectrum reduces to the 
ordinary Bethe-Heitler one. It is worth noting that in the diagrammatic language the 
above formula for N = 1 term corresponds to the set of diagrams shown in Fig. 1. We 
would like to emphasize that the spectrum (12), similarly to the gluon spectrum of Ref. 
[18], cannot be obtained if one ignores the logarithmic p-dependence of the function C(p) 
(10) at small p. One can see that, similarly to the gluon case, (11) and (12) do not contain 
any large logarithmic factor log(p.DPe//) (here p e // is the typical transverse scale) which 
one could expect there if one would neglect the p-dependence of C(p). 

The analytical formula (12) for the N = 1 spectrum demonstrates that due to the 
finite-size effects at large energies, when L~j > L, the induced spectrum is strongly peaked 
at x ~ 1. As will be seen below the N = 1 dominates for RHIC and LHC conditions, 
and the spectrum with the N > 2 scatterings included remains sharply peaked at x close 
to unity. Since the x-integral on the right-hand side of (1) is also dominated by large 
x one may expect that the induced photon emission to be an important source of the 
direct photons. It should be emphasized that our 1/(1 — x) spectrum does not mean 
enhancement of radiation at large x. It originates simply from suppression of photon 
emission at small x due to the finite-size effects. 

Now let us discuss the influence of the induced gluon emission on the photon 
bremsstrahlung ignored in the above derivation. The gluon emission ,i.e., processes like 
q — > 'jgq, and multi-gluon radiation, leads to quark energy loss, and should reduce radi- 
ation of hard photons. It is important that for soft gluons, dominating the quark energy 
loss, the gluon formation length, L 9 f ~ a 2 ^ f a"n _ > , turns out to be of the same order as 
that for photons. For this reason, as was already noted, photon and gluon emission should 
be treated on an even footing. An accurate analysis with gluon emission is a complicated 
task which is far beyond the scope of the present paper. We account for influence of the 
induced gluon emission at qualitative level. In line with the prescription suggested in [19] 
for hadron spectra, we replace in (1) the cross section of quark production by an effective 
cross section 

d <effM _ Nb f dx dI(x,p T /(l-x)) daf{p T /{l-x)) 
dydp\ m J (1 — x) 2 dx dydp\ 

where dI(x,E)/dx is the probability distribution in the quark energy loss, and Nu n is 
the number of the binary nucleon-nucleon collisions. In the p^-region of interest the cross 
section of quark production in pp collisions can be parametrized as 



dydpl ~ (pt + Po)" 
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withpo ~ 1-6 GeV, n ^ 8 for RHIC, andp ~ 0.6 GeV, n 5.3 for LHC, the normalization 
constant A will not be important to us. For such a ^-dependence to good accuracy one 
can write the effective cross section of quark production as 

d ^)jeff „ N jggfaj f A ^ _ ^ n(pT) _ 2 dI(x,p T ) 



dydp\ bm dydp\ 



fdx(i-xr^- 2 . 1[x ; PT) , (is) 

J dx 



where 
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Then the effective fragmentation function for photon radiation accounting for gluon emis- 
sion can be written in the form 



D e /Jjx,E) = S g (E) 



dP vac (x,E) dP ind (x,E) 



(16) 



dx dx 
where S g (E) is the gluon suppression factor given by 

S g (E) « P (E) + f dx(l - x )n(E)-2dI(^E) 

Here we separated the probability of photon emission without gluons, Po- In terms of the 
gluon spectrum, dP g /dx, it reads 

P (E) = exp (- I' dx dP ^ E ^ 



dx 



with x m i n ~ rrig/E. In numerical calculations we take x m i n = 2m g /E. 
We take the spectrum in the radiated energy entering (13) in the form 



dx \ Jx dy I dx 



The formula (17) is similar to the electron energy loss spectrum derived in [21]. It works 
well for small energy loss AE E. For RHIC and LHC the parametrization (17) 
reproduces the energy loss spectrum evaluated assuming independent gluon radiation 
with accuracy ~ 10 — 40%. It is enough to make qualitative estimate of the effect of gluon 
emission on the photon spectrum. An accurate calculation of the gluon suppression factor 
in the approximation of independent gluon emission [19, 22] does not make sense because 
this approximation itself does not have any serious theoretical justification. It should be 
noted that the suppression of hadron spectra evaluated using the energy loss spectrum 
(17) agrees well with that observed at RHIC. Thus, our approximation in some sense is 
justified by the experimental data. 

Neglecting the ISI effects the nuclear modification factor defined as 

td { \ 1 da^ A (p T )/dydp 2 T 

Raa{Pt) ~ W n ' daTM/dydpl (18) 
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can be approximately written as 

Jo Gferr™^- 2 [<iP mc (x,p T )/rfx + rfP ind (a;,p T )/cia;] 

**,<!*)« S.M /o w^w^, PT )/ fe ' (19 > 

We take the vacuum distribution in the form 

iP -f' E) = ^-(4 - 4, + 2S) ■ f™ dk'-JL-, , (20) 
dx Anx Jo (Ar + e z ) z 

where e = m q x, and A; maa , ps Pmax(x, (1 — x)). 

To evaluate the gluon suppression factor we use the formula for gluon spectrum 
similar to (3). However, we take the thickness equals L/2. It seems to be a reasonable 
choice to account for the fact that for L g f ~ about half of gluons are radiated after 
photon emission and cannot affect the photon bremsstrahlung. Note that this reduces the 
suppression effect of gluon emission as compared with the analyses [4, 5, 6, 7]. 

One remark is in order in connection with the above treatment of the gluon effects. 
Our gluon suppression factor includes only the induced gluon radiation. As far as the or- 
dinary vacuum hard gluon radiation is concerned, we assume that, due to small formation 
length, the corresponding suppression factors are approximately the same for AA- and 
pp-collisions. For this reason one may ignore the vacuum gluon radiation in evaluating 
the nuclear modification factor (if one uses the leading order vacuum spectrum (20)). 

3. For numerical calculations we use the one-loop running coupling constant with Aqcd = 

0. 3 GeV frozen at the value a s = 0.7. This parametrization is motivated by wanting to 
satisfy the relation 

r-2 GeV ry (h) 

dk ^± ~ o.36 GeV (21) 

obtained from the analysis of the heavy quark energy loss [23]. To fix the quark and Debye 
screening mass we use the results of the analysis within the quasiparticle model [24] of 
the lattice results. For the relevant range of temperature of the plasma T ~ (1 — 3)T C 
(T c 170 MeV is the temperature of the deconfinement phase transition) the analysis [24] 
gives for the quark and gluon quasiparticle masses m q ps 0.3 and m g ps 0.4 GeV. With the 
above value of m g from the perturbative relation = \f2m g one obtains /xp ~ 0.57 GeV. 
We assume the Bjorken [25] longitudinal expansion of the QGP with Tr 3 = T r 3 . For the 
initial conditions we use the values suggested in [20] To = 446 MeV and tq = 0.147 fm for 
RHIC, and T = 897 MeV and r = 0.073 fm for LHC. For RHIC the above condition were 
obtained from the charged particle pseudorapidity density dN/dy pa 1260 measured by 
the PHOBOS experiment [26] assuming an isentropic expansion and rapid thermolization 
at t ~ 1/3T . The LHC parameters correspond to dN/dy ~ 5625 estimated in [27]. 
Note that, since the dominating p-scale in (3) oc \fz for z <^ Lj, our results are not very 
sensitive to r . For the upper limit of the ^-integration in (21) we take L = Ra ~ 6 fm 

1 . This seems to be a reasonable value for central heavy-ion collisions since due to the 
transverse expansion the hot QCD matter should cool quickly at r > R A [25]. 

1 For our choice of the initial conditions the life-time of QGP is ~ 3 fm for RHIC. However, in the 
interval r <~ 3 — 6 fm the density of the mixed phase is practically the same as that for the pure QGP 
phase. 
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In Figs. 2 and 3 we show the x- dependence of the probability distribution of the 
induced q — > 75 transition (solid line) for several quark energies for RHIC and LHC. 
For comparison the vacuum spectrum (dashed line) is also shown. One sees that in the 
region of large x, which dominates the ^-integrals on the right-hand side of Eq. (19), the 
induced radiation exceeds the vacuum one (especially for the LHC case). To illustrate the 
influence of the LPM effect on the induced photon emission we also show the contribution 
of the N — 1 scattering (long-dashed line). It is seen that for x ~ 0.6 — 0.8 the LPM 
effect reduces the induced bremsstrahlung by a factor of ~ 0.8 for RHIC and ~ 0.5 for 
LHC. Note that the LPM effect and finite-size effects diminish in strength as x — > 1 since 
Lj becomes small and the spectrum should be close the Bethe-Heitler one in this limit. 
To demonstrate the role of the finite kinematic limits (neglected in (3)) in Figs. 2, 3 we 
also show the N = 1 scattering contribution evaluated using the set of diagrams shown in 
Fig. 1 with finite kinematic limits (dotted line). One can see that the kinematic effects 
become important at x close to unity. Numerical calculations show that they reduce the 
integral over x in the numerator in (19) by ~ 20 — 30%. However, this suppression to good 
accuracy is compensated by the increase in the gluon suppression factor due to similar 
kinematic effects for gluon emission. For this reason we neglect the kinematic effects in 
calculation of the nuclear modification factor (19). 

In Fig. 4 we plot results for the pr-dependence of the nuclear modification factor (19). 
We also show the results without gluon suppression factor (dashed line), i.e., for S g = 1. 
One sees that despite strong suppression due to gluon emission we obtain Raa ^ 1 
at pt ~ 5 — 15 GeV. This means that the FSI effects can enhance the direct photon 
production (especially for LHC energies, where bremsstrahlung is the dominating partonic 
mechanism). Of course, one should bear in mind that uncertainties in our theoretical 
predictions can be about 30-=-50%. Nevertheless, one can say that it is practically excluded 
that in central AA-collisions the direct photon production is strongly suppressed due to 
quark energy loss, as was predicted in [4, 5, 6, 7]. 
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Figures 





Figure 1 : The set of Feynman diagrams corresponding to the N — 1 scattering induced 
spectrum. 
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Figure 2: The photon spectrum of the q — > ^q transition for RHIC conditions. The solid 
line shows contribution from the induced photon emission calculated using the formula (3). 
The dashed line shows the vacuum spectrum (20). The N — l scattering contribution to 
the induced spectrum is shown by the long-dashed curves for infinite kinematic boundaries, 
and by the dotted curves for finite kinematic boundaries. 
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Figure 3: The same as in Fig. 2 but for LHC. 
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Figure 4: The p T -dependence of the nuclear modification factor (19) (solid line) for RHIC 
(a) and LHC (b) conditions. The dashed line shows the results without gluon suppression 
factor S g on the right-hand side of (19). 
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